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Biased Brownian motion of point-size particles in a three-dimensional tube with smoothly varying cross- 
section is investigated. In the fashion of our recent worki we employ an asymptotic analysis to the stationary 
probability density in a geometric parameter of the tube geometry. We demonstrate that the leading order 
term is equivalent to the Fick- Jacobs approximation. Expression for the higher order corrections to the 
probability density are derived. Using this expansion orders we obtain that in the diffusion dominated 
regime the average particle current equals the zeroth-order Fick- Jacobs result corrected by a factor including 
the corrugation of the tube geometry. In particular we demonstrate that this estimate is more accurate 
for extreme corrugated geometries compared to the common applied method using the spatially dependent 
diffusion coefficient D(x,f). The analytic findings are corroborated with the finite element calculation of a 
sinusoidal-shaped tube. 
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Particle transport in micro- and nanostructured 
channel structures exhibits peculiar characteris- 
tics which differs from other transport phenom- 
ena occurring for energetic systems. The the- 
oretical modelling involves Fokker-Planck type 
dynamics in three dimensions which cannot be 
solved for arbitrary boundary conditions imposed 
by the geometrical restrictions. Recently, much 
effort is drawn on a reduction of the complex- 
ity of the problem resulting in the so-called Fick- 
Jacobs approximation in which (infinitely) fast 
equilibration in certain spatial directions is as- 
sumed. Within the present manuscript we de- 
rive a reduction method which (i) corresponds in 
zeroth order in the expansion parameter, which 
describes the corrugation of the tube wall, to the 
celebrated Fick-Jacobs result and (ii) extends the 
validity of the Fick-Jacobs approximation towards 
extreme corrugated tube structures. 



I. INTRODUCTION 

The transport of large molecules and small parti- 
cles that are geometrically confined within zeolites 2 - - — , 
biological^ as well as designed nanoporcs£~— , channels 
or other quasi-one-dimcnsional systems attracted atten- 
tion in the last decade. This activity stems from the 
profitableness for shape and size selective catalysi o 9 ' 10 , 
particle separation and the dynamical characterization 
of polymers during their translocation^— . In partic- 
ular, the latter theme which aims at the experimental 
determination of the structural properties and the amino 
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acid sequence in DNA or RNA when they pass through 
narrow openings or the so-called bottlenecks, comprises 
challenges for technical developments of nanoscalcd chan- 
nel structures 1 ^ - — . 

Along with the progress of the experimental techniques 
the problem of particle transport through corrugated 
channel structures containing narrow openings and bot- 
tlenecks has give rise to recent theoretical activities to 
study diffusion dynamics occurring in such geometries 1 ^. 
Previous studies by Jacobs^ and Zwanzig 2 ^ ignited a 
revival of doing research in this topic. The so-called 
Fick-Jacobs approach^—, accounts for the elimination 
of transverse stochastic degrees of freedom by assum- 
ing a (infinitely) fast equilibration in those transverse 
directions. The theme found its application for biased 
particle transport through periodic 3D planar channel 
structure a 1 ' 21 " — as well as for tubesi2£~— exhibiting 
smoothly varying side- walls. 

Beyond the Fick-Jacobs (FJ) approach, which is suit- 
ably applied to channel geometries with smoothly vary- 
ing side-walls, there exist yet other methods for describ- 
ing the transport through varying channel structures like 
cylindrical septate channels 2 ^—, tubes formed by spher- 
ical compartment s 33 ' 34 or channels containing abrupt 
changes of cross diameter a 35 ' 36 . 

In a recent worki we have provided a systematic treat- 
ment by using a series expansion of the stationary prob- 
ability density in terms of a geometric parameter which 
specifies the channel corrugation for biased particle trans- 
port proceeding along a planar three-dimensional chan- 
nel exhibiting periodically varying, axis symmetric side- 
walls. We have demonstrated that the consideration of 
the higher order corrections to the stationary probability 
density leads to a substantial improvement of the com- 
monly employed Fick-Jacobs approach towards extreme 
corrugate channels. The object of this work is to pro- 
vide an analytic treatment to biased Brownian motion 
in cylindrical three-dimensional tubes with periodically 
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Figure 1. Sketch of a segment of a cylindrical tube with sinu- 
soidally varying radius R(x) that is confining the motion of 
the overdamped, point-like Brownian particle. The periodic- 
ity of the tube structures is L, the minimal and maximal tube 
widths are Aw and Afi, respectively. The constant force F 
pointing in the direction of the tube is applied on the parti- 
cles. 



varying radius. 

In Sec.|TI]we introduce the model system: a Brownian 
particle in a confined tube geometry with periodically 
modulated boundaries. The central findings, namely the 
analytic expressions for the probability density are pre- 
sented in Sec. IIIII It turns out that the latter scales 
linearly with the particle velocity derived within the FJ 
approach, cf. Sec. IIV1 In Sec. [V] we calculate the parti- 
cle mobility and the diffusion coefficient and employ our 
analytical results to a tube with sinusoidal varying cross- 
section. Section [VII summarizes our findings. 



II. TRANSPORT IN CONFINED STRUCTURES 

The present paper deals with biased transport of over- 
damped Brownian particles in a cylindrical tube with pe- 
riodically varying cross-section, respectively, radius R(x). 
A sketch of a tube segment with period L is shown in 
Fig. [1] The particles budge within a static fluid with 
constant friction coefficient r\. As the particle radius (as- 
sumed to be point-like) is small compared to the tube ra- 
dius, hydrodynamic particle-particle interactions as well 
as hydrodynamic particle-wall interactions can safely be 
neglected. Further the particles are subjected to an ex- 
ternal force with static magnitude F acting along the lon- 
gitudinal direction of the tube e^, i.e. the corresponding 
potential is U(x,r,(j)) = —Fx. 

The evolution of the probability density P (q, t) of 
finding the particle at the local position q = (x, r, </>) 
at time t is governed by the three-dimensional Smolu- 



fi t i>(q,f) + V q -J(q J t) = 0. 
where 



(la) 



J(q,t)= -P(q,t) e x -— V q P(q,i) (lb) 

T) 7] 

is the probability current J (q, t) = ( J x , J r , J^) T associ- 
ated to the probability density P (q, t) . The Boltzmann 
constant is ks and T refers to the environmental temper- 
ature. At the tube wall the probability current obeys the 
no-flux boundary condition (be) caused by the impene- 
trability of the tube walls, viz. J (q, t) ■ n = where n is 
the out-pointing normal vector at the tube walls. For a 
tube with radius R(x) the be read 



R'{x)J x (q, t) = J r (q, t) , for r = R{x) 



(2a) 



The prime denotes the derivative with respect to x. As 
a result of symmetry arguments the probability current 
must be parallel with the tube's centerline at r = 



J r (q,t)\ r 



= 0. 



(2b) 



Further the probability density satisfies the normal- 
ization condition J" it _ ri ^(l: *) <^ 3f l = 1 as well 
as the periodicity condition P(x + m L,r, (p, t) = 
P(x, r, (p, t) , Vto 6 Z. 

Since the external force acts only in longitudinal direc- 
tion P (q, t) is radial symmetric. This allows a reduction 
of the problem's dimensionality from 3D to 2D by inte- 
grating Eq. (jla[) over the angle </>; yielding 



dtP{x,r,t) 



d x 



k B T 1 



T d x (e «b t P(x,r,t) 
d r [rd r P{x,r,t)] , 



rj r 

where the two-point probability density is defined as 

2tt 

1 



(3) 



P(x,r,t)=— I d(t>P(x,r,frt) 



(4) 



Integrating Eq. ([3]) further over the cross-section and 
taking the boundary conditions Eqs. ([2]) into account, 
one gets 



d t P(x,t) J^d x 



R(x) 



dr r 



e "b t d x (e "b t P(x,r,t)j 
(5) 



Thereby the marginal probability density reads 

R(x) 2tt 



P(x,t) 



1 

= 2^ 



dr ', 



>P(x,r,4>,t) . 



(G) 
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In Rcfi we present a perturbation series expansion in 
terms of a geometric parameter for the problem of biased 
Brownian dynamics in a planar three-dimensional chan- 
nel geometry. Below we apply this method for Brownian 
motion in cylindrical three-dimensional tubes. In doing 
so we introduce dimensionless variables. We measure the 
longitudinal length in units of the period length L, viz. 
x = x/L. For the rescaling of the r-coordinate, we intro- 
duce the dimensionless aspect parameter e, i.e. the dif- 
ference of the widest cross-section of the tube, i.e. AO, 
and the most narrow constriction at the bottleneck, i.e. 
Aw, in units of the period length, yielding 



(Afi - Aw) 



At steady state, the Smoluchowski equation Eq. (|3|) in 
dimensionless units becomes 

e 2 d x [e- u d x (e u p st {x,r))] + \ [rd rPst (x,r)) = 0, 

(11) 

and the no-flux boundary conditions Eqs. ([2]) read 

= [d r p st (x,r) - e 2 h'(x)e- u d x (e u p st (x,r))] r=h(x) , 

(12a) 

=d r p st (x, r)| r=0 . (12b) 



(7) 



The dimensionless parameter e characterizes the devi- 
ation of the boundary from the straight tube corre- 
sponding to e = 0. Several authors considered different 
choices for the expansion parameter like the averaged half 
widt h 25 ' 39 ! 40 or the ratio of an imposed anisotropy of the 
diffusion constants e 2 = Dy/Dj^. 

We next measure, for the case of finite corrugation 
e 7^ 0, the radius r in units of eL, i.e. r = eLr 
and, likewise, the boundary function R(x) = eLh{x). 
Time is measured in units of r = L?r\j '(fc^T) which is 
twice the time the particle requires to overcome diffu- 
sively, at zero bias F = 0, the distance L, i.e. t = t/r. 
The potential energy is rescalcd by the thermal en- 
ergy k^T, i.e., for the considered situation with a con- 
stant force component in longitudinal direction, one gets 
U := U(x, r) = —Fx/{k&T) = — fx, with the dimension- 
less force magnitud o 21 ' 23 : 



/ = 



FL 



(8) 



After scaling the probability distribution reads p (q, t) = 
e 2 L 3 P (q, t) . Below we shall omit the overbar in our 
notation. 

Further we concentrate only on the steady state, i.e., 
limt_>. 00 p(x, r, t) := p at (x,r), which is in fact, the only 
state necessary for deriving the key quantities of particle 
transport like the average particle velocity (x) 



III. ASYMPTOTIC ANALYSIS 

We apply the asymptotic analysis^ - to the problem 
stated by Eq. (TlTj) and Eqs. p^|) . In doing so, we use for 
the stationary probability density p st (x,r) (the index st 
will be omitted in the following) the ansatz 



P{x,r) = ^2e 2n p n (x,r) 



(13) 



n=0 



in the form of a formal perturbation series in even orders 
of the parameter e. Substituting these expressions into 
Eq. <jTTJ), we find 



1 °° f 1 

= -d r [rd rPo (x, r)} + £2n \ ~ d r [rd rPn (x, r)\ 

n=l 

+d x [e- u d x (e u p n _ 1 (x,r))]} . 



(14) 



The no-flux be at the tube walls r = h(x), cf. Eq. (|12al) . 
turns into 



= d r po{x, r)+ Y £2n {d r Pn(x, r) 



(15a) 

-h'(x)e- u d x (e u p n . 1 (x,r))} , 

and the be at the center line of the tube r = 0, cf. 
Eq. (|12b|) . then reads 



= lim 



(x(t)) 



t— too t 



1 h(x) 

dx J drrJ*(x,r), (9) 
o o 



and the effective diffusion coefficient D c g in force direc- 
tion. The latter is given by 



-D e ff = lim 

t— J-oo 



(x 2 (t)) - (x(t)Y 
2t 



(10) 



and can be calculated by means of the stationary proba- 
bility density p st (x, r) using an established method taken 
from Refi 4 ^. 



0=J2z 2n drPn(x,r) 



n=0 



(15b) 



We claim that the normalization condition for the prob- 
ability density p(x, r) corresponds to the zeroth solution 
Po{x,r) that is normalized to unity, 

1 h(x) 

(p (x,r)) = J dx J drrp (x,r) = 1 . (16) 
o o 

Consequently the higher orders in the perturbation se- 
ries have zero average, (p n (x,r)) = 0,n = 1,2,3,... . 
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Further each order p n has to obey the periodic boundary 
condition p n (x + m, r) = p n (x, r) , Vm € Z. 

InSec. IIII A| we demonstrate that the zeroth order 
of the perturbation series expansion coincides with the 
Fick-Jacobs equatio n 19 ' 20 . Referring to 21 ' 42 an expres- 
sion for the average velocity (x} is known. Moreover, in 
Sec. IIII FJl the higher order corrections to the probability 
density are derived. Using those results we are able to 
obtain corrections, see in Sec. IV Bl to the average veloc- 
ity beyond the zeroth order Fick-Jacobs approximation 
presented in the next section. 



A. Zeroth Order: the Fick-Jacobs equation 



For the zeroth order, Eq. (jl4j) read 



1 



= -d r [rd r po{x,r)} 



(17a) 



supplemented with the corresponding no-flux boundary 
condition at r = as well as at r = h(x) 



=d r p (x,r) . 



(17b) 



We make the ansatz po(x,r) = g(x)e~ u where g(x) is 
an unknown function which has to be determined from 
the second order O (e 2 ) balance given by Eq. (|14|) : 

=d x (e- u g'(x)) + U r [rd r (pi(x.r))] . (18) 

Integrating the latter over the radius r and taking the 
no-flux boundary conditions Eqs. ([15} into account, one 
immediately obtains 

0=d x (e- A ^g'(x)) . (19) 

The effective entropic potential A(x) is defined by 



h ( x 

-A(x) = J 



drre 



-U(x,r) 



(20) 



and for the problem at hand the latter looks explicitly 

'h 2 (x) 



A(x) = — f x — In 



21) 



Note, that upon an irrelevant additive constant, i.e. 
ln(27r), the effective entropic potential corresponds to 
that given in Ref.— . 

The normalized stationary probability density within 
the zeroth order reads 3 -! 



e -U(x,r) X j e A(x') dx l 



Po(x,r) = e u g{x) =- 



x+l 

/ dxe-^x) J e A(x')d x i 

x 



, (22) 



and, moreover, the marginal probability density Eq. © 
becomes 



Po{x) 



-A(a 



9i x ) 



(23) 



Expressing next g(x) by po(x), see Eq. (|19p . then yields 
the celebrated stationary Fick-Jacobs equation 

} d x (e A ^p (x))} (24) 



= d x 



-A(x 



derived previously in Re h 20 ' 43 . 

Thus, we demonstrate that the leading order term of 
the asymptotic analysis is equivalent to the FJ-equation. 
In the FJ equation the problem of biased Brownian dy- 
namics in a confined 3D geometry is replaced by Brow- 
nian motion in the tilted periodic one-dimensional po- 
tential A(x). In general, the stationary probability den- 
sity of finding an overdamped Brownian particle budging 
in a tube with periodically varying cross-section is suf- 
ficiently described by Eq. (f2"4")l as long as the extension 
of the bulges of the tube structures is small compared to 
the periodicity, i.e. e< 1. 

Then the average particle current is calculated by in- 
tegrating the probability flux Jq over the unit-cel l 42 ' 44 ' 45 

1 h(x) 

W))o = dx drrJ$(x,r) 



o o 



1 - e~f 



J dxe A ( x ) J e- A ( x ')dx' 

x-l 



(25) 



In the spirit of linear response theory, the mobility in 
dimensionless units is defined by the ratio of the mean 
particle current Eq. (|25p and the applied force /, yielding 



Mo (/) = 



(Hf))o 
f 



(26) 



Note, that in order to obtain the mobility in physical 
units one has to multiply with the mobility of uncon- 
fincd particles, i.e. l/rj. 

Further, resulting from the normalization condition 
Eq. (Til)]) , the average particle velocity Eq. ([9]) simplifies 
to 



(x) = (x n - 



(d x p n {x,r)) 



(27) 



Therefore we derive that the average particle current is 
composed of (i) the Fick-Jacobs result (x) Q , cf. Eq. (|2"5j) . 
and (ii) becomes corrected by the sum of the averaged 
derivatives of the higher orders p n {x,r). We next ad- 
dress the higher order corrections p n (x,r) of the prob- 
ability density which become necessary for more corru- 
gated structures. 
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B. Higher order contributions to the Fick-Jacobs equation 

According to Eq. (fl"4)) , one needs to iteratively solve 

-d r [rd r p n (x,r)] =£p n -!(x,r) , n>l, (28) 
r 

under consideration of the boundary conditions Eqs. (|15l) . 
In Eq. (|28p . we make use of the operator £, reading 
£ = (fd x — d x ). Each solution of the second order 
partial differential equation Eq. (f28f possesses two in- 
tegration constants d„_i and d n> 2- The first one, d n _i, 
is determined by the no-flux be at the centerline r = 0, 
cf . Eq. (|15b[) , while the second provides the zero average 
condition (p n (x,r)) = 0, n > 1. 

For the first order correction, the determining equation 
reads 

U r [rd rPl {x,r)\ = 2(x) d x (r^j , (29) 
and after integrating twice over r, we obtain 

Pl (x,r) = - (x) Q (^v*. (30) 

Hereby as requested above, the first integration con- 
stant di,i(x) is set to in order to fulfill the no-flux 
bes, and the second must provide the normalization con- 
dition Eq. HH), i.e. d\ t 2 = 0. One notices that the 
first correction to the probability density becomes posi- 
tive if the confinement is constricting, i.e. for h'(x) < 
and (x) Q 7^ 0. In contrast, the probability density be- 
comes less in unbolting regions of the confinement, i.e. 
for h'(x) > 0. Please note, that the first order cor- 
rection scales linearly with the average particle current 
(x) Q . Overall, the break of spatial symmetry observed 
within numerical simulations in previous work o 22 ' 46 is re- 
produced by this very first order correction. Particularly, 
with increasing forcing, the probability for finding a par- 
ticle close to the constricting part of the confinement in- 
creases, cf. Refi^i 4 ^. 

Upon recursively solving, the higher order corrections, 
n > 1, are given by 

p n (x,r) =£ n p (x,r) — — + d„ >2 , (31) 
(2 n n\) 

Jin / i \ 

(2"n!) 2 \h 2 (x)J 

where the operator £ applied n-times yields the expres- 
sion 

" / \ pin+k 

r = E I Hfr^. (32) 

k=0 V / 

Note that each single order p n (x, r), cf. Eq. (f3"Tj) . satisfies 
the normalization condition, the be at the centerline but 



does not obey the be at the tube wall, see Eq. (|15a[) . 
The stationary probability density p(x, r) is obtained by 
summing all correction terms, cf. Eq. (fT3j) . yielding 

p(x,r) =p (x,r) + £ n p (x,r)- -j + d„, 2 ■ 

(33) 

Inserting Eq. ()33p into the equation for the no-flux be at 
the tube wall, cf. Eq. (|12ap . results to 

h(x) 

= Y. £2n+2dx / drre- u{x ^d x (e u ^ x ^ ) p n {x,r) S j 

h(x) 

=e 2 d x J drre- u{x ^ ) d x (e u{x ' r p{x 1 r) s j . (34) 
o 

According to Eq. © the latter equals zero in the steady 
state. 

Summing up, the exact solution for the stationary 
probability density of finding a biased Brownian particle 
in tube is given by Eq. (|3"3"|) . The latter solves the cor- 
responding Smoluchowski equation Eq. under satis- 
faction of the normalization as well as the periodicity re- 
quirements. More importantly the solution, cf. Eq. (|33|) . 
obeys the no-flux boundary conditions at the centerline 
r = as well as at the tube wall r = h(x). Further one 
notices that p(x, r) is fully determined by the Fick-Jacobs 
results po (x, r). Caused by £po(x,r) cx (i(/)) , the con- 
tribution of the higher order corrections to the 2D prob- 
ability density scales linear with the average particle cur- 
rent in the FJ limit (i) , cf. Eq. (f3"Tj) . The latter is deter- 
mined by the break of spatial symmetry induced by the 
external force. Consequently, in the absence of the exter- 
nal force / = the stationary probability density equals 
the zeroth order contribution p(x,r) = po(x,r) = const 
despite the value of e. 

According to Eq. (|27|) , one recognizes that the average 
particle current scales with the average particle current 
obtained from the Fick-Jacobs formalism (x) for all val- 
ues of e. Therefore, in order to validate the obtained 
results for p{x, r) as well as to derive correction to the 
mean particle current it is required to calculate (x) Q first. 

IV. TRANSPORT QUANTITIES FOR A SINUSOIDALLY 
SHAPED TUBE 

In the following we study the key transport quantities 
like the particle mobility /i(/) and the effective diffusion 
coefficient D e ff(/) of point-like Brownian particles mov- 
ing in a sinusoidally-shaped2£ tube. The dimensionless 
boundary function h(x) reads 

h(x)=-[-—:+sm(2Trx)), (35) 
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and is illustrated in Fig. [TJ The function h(x) is solely 
governed by the aspect ratio of the minimal and maxi- 
mum tube width 6 = Auj/AH. Obviously different re- 
alizations of tube geometries can possess the same value 
of S. The number of orders have to taken into account 
in the perturbation series Eq. (|13[) . respectively, the ap- 
plicability of the Fick-Jacobs approach to the problem, 
depends only on the value of the geometric parameter 
e = Ail (1 — 5) /L for a given aspect ratio 5. 

First, we obtain the particle mobility fio within the 
zcroth-order (Fick-Jacobs approximation). Referring to 
Eqs. (f2l)j) and ([25)1 the dimensionless mobility is given by 



1 - 

Mo(/J = — i 

/ / dxe~f x ' /h 2 (x>) J ef x h 2 {x)dx' 

x-l 



(36) 



For the considered sinusoidal boundary function, 
cf. Eq. (|35|) . we obtain 



1 



1 



Mo(/) 2 V / F _ T 



b(2b 2 



1 - 



46 P 



2VP~T -2b 3 + 3b) f 



P + (2nY 

(37) 



P + (^) 

with the substitute b = (1 + S)/(l — 5). Caused by the 
reflection symmetry of the boundary function h(x) the 
particle mobility obeys Ho(—f) = Mo(/)- Thus it is suffi- 
cient to discuss only the behavior for / > 0. 

In the limiting case of infinite large force strength, the 
mobility goes to 



/ 



lim fj, (f) =1 . 



(38) 



With decreasing force magnitude / the mobility de- 
creases as well till /io(/) attains the asymptotic value 



2VS 86 



(39) 



In the diffusion dominated regime, |/| <C 1, the 
Sutherland-Einstein relation emerge n 47 ' 48 and thus the 
dimensionless mobility equals the dimensionless effective 
diffusion coefficient: 



lim iu(f) = lim iD eff (f). 



(40) 



In the limit of vanishing bottleneck width, i.e. 6 — > 0, 
the mobility, respectively, D e g tends to 0. In contrast, 
for straight tubes corresponding to S = 1, i.e. e = 0, the 
mobility as well as the effective diffusion coefficient equal 
their free values which are one in the considered scaling. 

In a recent work 1 -, we studied the biased Brownian mo- 
tion in a planar three-dimensional channel geometry with 
periodically varying width 2h(x). We found that the 
asymptotic value for the mobility is given by 2^/(1 + 5) 
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Figure 2. (Color online) The particle mobility (a) and the 
effective diffusion constant (b) for a Brownian particle mov- 
ing inside a sinusoidal tube are depicted as function of the 
external force magnitude /. The maximum tube width is 
kept fixed, viz. Afi = 0.1, while the aspect ratio is varied 
8 — 0.01, 0.1, 1, respectively, the corresponding values for e 
are e = 0.099, 0.09, 0. The symbols correspond to the nu- 
merical obtained mobility, respectively, the effective diffusion 
coefficient. In subfigure (a) the lines correspond to the ana- 
lytic result, cf. Eq. PTl) . 



for / — > 0, cf. Eq. (45) in RefA Comparing this result 
with the asymptotic value Eq. ([39]) one notices that the 
mobility and the effective diffusion coefficient in a tube, 
respectively, are less compared to case of a planar channel 
geometry. 

In Fig. [5] we depict the dependence of the mobility 
/i(f) and the effective diffusion coefficient D e s(f) on the 
external force magnitude /. The numerical results are 
obtained by solving the stationary Smoluchowski equa- 
tion Eq. (fTaj) using finite element method 4 ^ and subse- 
quently calculating the average particle current accord- 
ing to Eq. ([9]). In order to determine the effective diffu- 
sion coefficient D e s(f), one has to solve numerically the 
reaction-diffusion equation for the B-fiel d 39 ^ 41 . 

Referring to Fig. [Ha) on e notices that the analytic 
predictions for the particle mobility Eq. (|37p arc corrob- 
orated by numerics. Further, one observes that for the 
case of smoothly varying tube geometry, i.e. Ail/L <C 1, 
the analytic result is in excellent agreement with the nu- 
merics for a large range of dimensionless force magni- 
tudes /, indicating the applicability of the Fick-Jacobs 
approach. As long as the extension of the bulges of the 
tube structures is small compared to the periodicity, suf- 
ficiently fast transversal equilibration, which serves as 
fundamental ingredient for the validity of the Fick-Jacobs 
approximation is taking place. 

The effective diffusion coefficient D c g(f) exhibits 
a non-monotonic dependence versus the dimensionless 
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Figure 3. (Color online) The influence of the geometric pa- 
rameter e on the particle mobility is presented. The value of 
e is varied, viz. e — 0.05,0.5,2.5 , while the aspect ratio is 
kept fixed, viz. S = 0.5. The solid line corresponds to the 
analytic result Eq. (I37|l while the dash-dotted line indicates 
the asymptotic value one. 



force /, see Fig. &h). It starts out with a value that is 
less than the free diffusion constant in the diffusion dom- 
inated regime, i.e. |/| <C 1. According to the Suthcrland- 
Einstcin-rclation the value equals the mobility value, 
cf. Eq. (|39p . Then it reaches a maximum with increasing 
/ and finally approaches the value of the free diffusion 
from above. Further one notices that the location of the 
diffusion peak as well as the peak height depends on the 
aspect ratio S. With decreasing width at the bottleneck, 
while keeping the maximum width AO fix, the diffusion 
peak is shifted towards larger force magnitude /. Simul- 
taneously the peak height grows. In the limit of a straight 
tube, i.e. 5 — > 1, as expected the effective diffusion co- 
efficient coincide with its free value which is one in the 
considered scaling. 



In Fig. [3] we present the impact of the expansion pa- 
rameter e on the particle mobility. It turns out that for 
values of e < 0.1 the Fick- Jacobs approximation is in very 
good agreement with the simulation. With increasing ge- 
ometric parameter e the difference between the FJ-result 
(solid line) and the numerics is growing. The more avail- 
able space in the tube leads to a decrease of the particle 
mobility, respectively, of the effective diffusion coefficient. 
Consequently, the higher order corrections to the station- 
ary probability density p{x, r), see Eq. (|13|) . respectively, 
the corrections to the mobility, cf. Eq. (j2"T)) . need to be 
included in order to provide a better agreement. 



V. CORRECTIONS TO THE MOBILITY AND 
DIFFUSION COEFFICIENT 

A common used way to include the corrugation of the 
channel structure bases on the concept of the spatially 
dependent diffusion coefficient D(x, /) which was intro- 
duced by Zwanzig2£ and subsequently supported by the 
study of Reguera and Rubi 4 ^. Zwanzig obtained the FJ 
equation, cf. Eq. (|2"4")) . from the full 3D Smoluchowski 
equation upon eliminating the transverse degrees of free- 
dom supposing infinitely fast relaxation. In a more de- 
tailed view, we have to notice that diffusing particles can 
flow out from/ or towards the wall in r-direction only 
at finite time. These finite relaxation processes are in- 
cluded by scaling the diffusion constant in longitudinal 
direction by a position dependent function, viz. D(x, /). 
The latter substitutes the constant diffusion coefficient - 
which is one in the considered scaling - in the common 
stationary FJ equation, cf. Eq. ([24]) , yielding 



= flL 



D(xJ)e- A tod x (e A Wp(3 



(41) 



According to Eq. (|4"Tj) an expression for the particle mo- 
bility that is similar to the Stratonovich formul a 21 ' 42 for 
the mobility in titled periodic energy landscapes, but now 
including the the spatial diffusion coefficient D(x, /), can 
be derive d 21 ' 22 . In the diffusion dominated regime, i.e. 
| /| <C 1, this expression simplifies to the Lifson- Jackson 
formul a 22 ! 50 



lim fx(f) = lim D cS (f) = 

(h 2 (x)) 



D(xfi)h' 2 {x) 



(42) 



with the period average (•) = f Q -dx. Unfortunately, for 
many boundary functions h(x) it is impossible to analyt- 
ically evaluate the expression Eq. (|4"2"]l . 



A. Spatially dependent diffusion coefficient D(x, f) 

A first systematically treatment taking the finite dif- 
fusion time into account was presented by Kalinay and 
Percus (KP ) 24 i 51 . Their suggested mapping procedure 
enables the derivation of higher order corrections in terms 
of an expansion parameter s 2 KP , which is the ratio of the 
diffusion constants in the longitudinal and transverse di- 
rections. Within this scaling, KP have shown that the 
fast transverse modes (transients) separate from the slow 
longitudinal ones and therefore the transients can be pro- 
jected out by integration over the transverse directions. 

In what follows we present a derivation for the spatially 
dependent diffusion coefficient D(x, f) which based on 
our previously considered perturbation series expansion 
for the stationary probability density Sec. IIIII Accord- 
ing to Eq. ([5]) the marginal probability current J x (x), 
equivalent to Eq. (|41[) , can be derived in an alternative 
way using the stationary two-point probability density 
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p(x,r), yielding 

-J x {x) =D(x,f)e- A ^d x (e A ^p{x) 



h{x) (43) 

drre- u{x ' r) d x (e u ^p(x,r)j . 



The second equality determines the sought-after spatial 
dependent diffusion coefficient D(x,f). 

One immediately notices that the relation Eq. (|43|) sim- 
plifies to D(x, f) f p(x) = dr r f p(x, r) in the force 
dominated regime |/| ^> 1. Then it follows that the spa- 
tially diffusion coefficient equals the free one, which is 
one in the considered scaling, 



lim D(x,f) = 1 

f-t-oo 



(44) 



In the opposite limit of small force strengths, i.e. for 
|/| <C 1, diffusion is the dominating process. Then 
Eq. simplifies to 



D(xJ) h 2 (x)d x 



h 2 {x) 



h(x) 



dr r d x p (x, r) . (45) 



Inserting our result for the stationary probability density 
p(x,r), cf. Eq. (|33j) . into Eq. (|45|) we determine an ex- 
pression for D(x, /). In compliance with Ref^, we make 
the ansatz that all but the first derivative of the boundary 
function h{x) are negligible. Then the n-times applied 
operator £, cf. Eq. ((32]), simplifies to £" = (-1)" d "~ 
yielding, 

p n (x, r) = 2 (-1)" (x) £ ";; 2 [ll J n+1 r 2 " + 0{h"{x)) 



9i 2 » ' 



(2"n!) 2 h 2n+1 



(46) 



Inserting the latter into Eq. (|43|) and calculating the com- 
plete sum, one finds 



lim D(x, f) = — ^^^=^^^= 

yjl + (eh'(x)f 



+ 0(h"(x)) (47) 



for the spatially dependent diffusion coefficient D(x, f) in 
the diffusion dominated regime, i.e. for |/| <C 1. Using 
our above presented series expansion for the stationary 
probability density p{x, r) we confirm the expression for 
the spatially dependent diffusion coefficient D{x, f) pre- 
viously derived by Reguera and Rubi^ and K P 24 i 51 . 



B. Corrections based on perturbation series expansion 

Next, we derive an estimate for the mean particle cur- 
rent (x(f)) based on the higher expansion orders p n {x, r). 
Referring to Eq. (1271) , the average particle current is com- 
posed of (i) the Fick- Jacobs result (i) , cf. Eq. (f2"5")l . and 



(ii) becomes corrected by the sum of the averaged deriva- 
tives of the higher orders p n (x,r). Immediately one no- 
tices that the integration constant d n _2, resulting from 
the normalization condition Eq. (|I6p. does not influence 
the result for the average particle velocity, cf. Eq. (f27|) . 

We concentrate on the diffusion dominated limit 
|/| I and further we make the ansatz that all but 
the first derivative of the boundary function h(x) are 
negligibl o 1 ' 24 . Then the partial derivative of p n (x, r) with 
respect to x simplifies to 

d xPn (x, r) = 2 (x) (-r ^^y + 0(h"(x)) 



(2 n n\y 



(48) 



Inserting the latter into Eq. (|2"T|) and integrating the re- 
sults over one unit-cell, results in 

(- 1 ) n ( 2 -+ 1 )' V(,)) 2 " 



71 = 



(n!) (n + 1)! 



yJl + (eh'(x)Y 



(49) 



We obtain that the average transport velocity is obtained 
as the product of the zeroth-ordcr Fick- Jacobs result and 
the expectation value of a complicated function includ- 
ing the slope of the boundary for |/| <C I. In the pre- 
viously studied case of biased Brownian motion in a 3D 
planar channel geometry^ we have found that the aver- 
age transport velocity is obtained as the product of the 
zeroth-ordcr Fick- Jacobs result and the expectation value 
of the spatially dependent diffusion coefficient (D(x,0)). 
In contrast to the 3D planar geometry, for biased trans- 
port in extreme corrugated tubes the corrections term to 
the particle velocity does not coincide with the expecta- 
tion value of D(x,0), see Eq. (T47]) . 

Calculating the expectation value in Eq. (j49]l for the 
considered tube geometry, cf. Eq. (|3"5"|) . yields to the 
estimate 



lim <*(/)>* lim( i(/) ) ^(i,i, 2, -(f 



(50) 



where %Fi (•) is the first hypergeometric function. We 
derive that in the diffusion dominated regime the aver- 
age velocity is obtained as the product of the zeroth-order 
Fick- Jacobs result and the correction term 2F1 (•) includ- 
ing the corrugation of the tube structure. Referring to 
the Sutherland-Einstein relation, cf. Eq. (|40|) . if the av- 
erage current (x(f)) (or the effective diffusion coefficient 
£>° ff (/)) is known in the zeroth order, the higher order 
corrections to both quantities can be obtained according 
to Eq. ([5fJ|). 

We have to emphasize that considering only the 
first derivative of the boundary function h'(x) 
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Figure 4. (Color online) Comparison of the analytic the- 
ory versus precise numerics (in dimensionless units): The 
mobility and the effective diffusion constant for a Brownian 
particle moving inside a tube with sinusoidal varying radius 
are depicted as function of geometric parameter e in units 
of the maximum channel width AQ. The latter is varied 
Afl/L = 0.1,1,5 (from top to bottom) while the external 
bias is kept fixed / = 10 -3 (corresponding to the diffusion 
dominated regime). The symbols correspond to the numerical 
obtained mobility (triangles) and the effective diffusion coef- 
ficient (circles). The solid lines correspond to analytic higher 
order result, cf. Eq. (|50p . The zeroth order - Fick- Jacobs re- 
sults given by Eq. ((39} collapse to a single curve hidden by 
the solid line for AQ./L = 0.1. In addition the numerical 
evaluation of Eq. (I42|) is represented by the dash-dotted lines. 



results in an additional term proportional to 

e 2 2-Fi (3/2,3/2,3,- (£7r/2) 2 ). Taking further the 

second derivative h"(x) into accounts indicates that this 
second term is negligible compared to 2F1 (1/2, ■ ■ •) for 
arbitrarily value of e. 

In Fig. 31 we present the dependence of the /Lt(/) (tri- 
angles) and D e g(f) (circles) on the slope parameter e for 
/ = 10~ 3 . One observes that the numerical results for 
the effective diffusion coefficient D c g(f) and the mobil- 
ity coincide for all values of e, thus corroborating 
the Sutherland-Einstein relation. In addition, the Fick- 
Jacobs result, given by Eq. (|39[) . the higher order result 
(solid lines), see Eq. ([50)) . and the numerical evaluation 
of the Lifson- Jackson formula using D(x, 0) (dash-dotted 
lines), cf. Eq. (T4"2")) . arc depicted in Fig. HI 

For the case of smoothly varying tube geometry, i.e. 
AQ/L <C 1, all analytic expressions are in excellent 
agreement with the numerics, indicating the applicabil- 
ity of the Fick- Jacobs approach. In virtue of Eq. ([7]) , the 
geometric parameter is defined by e = (A£l — Aoj) /L 
and hence the maximal value of e equals ACt/L. Con- 
sequently the influence of the higher expansion orders 
e 2n (d x p n (x,r)) on the average velocity Eq. (|2"T|) and on 
the mobility, respectively, becomes negligible if the max- 



imum tube's width AQ/L is small. 

With increasing maximum width the difference be- 
tween the F J-result Eq. (|3"5j) and the numerics is growing. 
Specifically, the FJ-approximation overestimates the mo- 
bility p and the effective diffusion coefficient -D ff- Con- 
sequently the corrugation of the tube geometry needs to 
be included. The consideration of D(x,Q), cf. Eq. (f4"2"l) 
provides a good agreement for a wide range of e-values 
as long as the maximum width AQ/L is on the scale to 
the period length of the tube, i.e. AVL/L ~ 1. Upon 
further increasing the maximum width AD, / L diminishes 
the range of applicability of the presented concept. In 
detail the expression Eq. (f42j) drastically underestimate 
the numerical results due to the neglect of the higher 
derivatives of the boundary function h(x). Put differ- 
ently, the higher derivatives of h(x) become significant 
for AQ/L > 1. 

In contrast, one notices that the result for the mobil- 
ity and the diffusion coefficient, respectively, based on 
the higher order corrections to the stationary probability 
density Eq. (|50[) is in very good agreement with the nu- 
merics. For tube geometries where the maximum width 
AQ./L is on the scale to the period length, i.e. AVl/L ~ 1, 
the correction estimate matches perfectly with the nu- 
merical results. Further increase of the tube width results 
in a small deviation from the simulation results. 



VI. SUMMARY AND CONCLUSION 

In summary, we have considered the transport of point- 
sized Brownian particles under the action of a constant 
and uniform force field through a 3D tube. The cross- 
section, respectively, the radius of the tube varies peri- 
odically. 

We have presented a systematic treatment of particle 
transport by using a perturbation series expansion of the 
stationary probability density in terms of a smallness pa- 
rameter which specifies the corrugation of the tube walls. 
In particular, it turns out that the leading order term of 
the series expansion is equivalent to the well-established 
Fick-Jacobs approac h 19 1 20 . The higher order corrections 
to the probability density become significant for extreme 
bending of the tube's side- walls. Analytic results for each 
order of the perturbation series have been derived. Sim- 
ilar to biased Brownian motion in a 3D planar channel 
all higher order corrections to the stationary probability 
density scale with the average particle current obtained 
from the Fick-Jacobs formalism. 

Moreover, for the diffusion dominated regime, i.e. 
for small forcing |/| <C 1, we calculate the correction 
to the mean particle velocity originated by the tube's 
corrugation using the series expansion for the station- 
ary probability density. According to the Sutherland- 
Einstein relation, the obtained relation is also valid for 
the effective diffusion coefficient. In addition, by using 
the higher order corrections we present an alternative 
derivation for the spatially dependent diffusion coefficient 
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D(x, /) which substitutes the constant diffusion coeffi- 
cient present in the common Fick- Jacobs equation based 
on similar assumptions as those suggested by Kalinay and 
Percus as well as by Rubi and Reguera. 

Finally, we have applied our analytic results to a spe- 
cific example, namely, the particle transport through a 
tube with sinusoidally varying radius R(x). We cor- 
roborate our theoretical predictions for the mobility and 
the effective diffusion coefficient with precise numerical 
results of a finite element calculation of the stationary 
Smoluchowski-cquation. 

In conclusion, the consideration of the higher order 
corrections leads to a substantial improvement of the 
Fick- Jacobs-approach, which corresponds to the zeroth 
order in our perturbation analysis, towards more wind- 
ing boundaries of the tube. Notably, we have shown that 
the common approach using the the spatially dependent 
diffusion coefficient D(x,f) fails for extreme corrugated 
tube geometries. 
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